Noname manuscript No. 

(will be inserted by the editor) 



Ball on a Beam: Stabilization under Saturated Input 
Control with Large Basin of Attraction 

Yannick Aoustin* and Alexander Formal'skii 



the date of receipt and acceptance should be inserted later 



Abstract This article is devoted to the stabilization of two under actuated planar 
systems, the well known straight beam-and-ball system and an original circular beam- 
and-ball system. The feedback control for each system is designed, using the Jordan 
form of its model, linearized near the unstable equilibrium. The limits on the voltage, 
fed to the motor, are taken into account explicitly. The straight beam-and-ball system 
has one unstable mode in the motion near the equilibrium point. The proposed control 
law ensures that the basin of attraction coincides with the controllability domain. The 
circular beam-and-ball system has two unstable modes near the equilibrium point. 
Therefore this device, never considered in the past, is much more difficult to control 
than the straight beam-and-ball system. The main contribution is to propose a simple 
new control law, which ensures, by adjusting its gain parameters, that the basin of 
attraction arbitrarily can approach the controllability domain for the linear case. For 
both nonlinear systems, simulation results are presented to illustrate the efficiency of 
the designed nonlinear control laws and to determine the basin of attraction. 
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1 Introduction 

Among the mechanical systems, the under actuated systems, which have fewer con- 
trols than configuration variables, represent a great challenge for the control. They are 
characterized by the underactuation degree, which is the difference between the num- 
bers of configuration variables and controls. An active field of research exists, due to 
the applications of under actuated systems such as aircrafts, satellites with thrusters, 
spacecrafts, flexible robots, legged robots, which adopt a dynamical stable walking 
or running gait, inverted pendulums. For example, for a planar vertical take-off and 
landing aircraft (PVTOL), an approximate input-output linearization procedure is de- 
veloped in pP to get a bounded tracking and an asymptotic stability. In paper [2], the 
stabilization of a satellite is studied when one of its three thrusters is not efficient. 
Flexible robots have an infinite number of flexible modes, which can be damped, using 
a controller, based on a discrete model, (see [3], [4], etc.). Inverted pendulum devices 
are used like a testbed research or for education application to investigate new con- 
trol laws for the stabilization or the swing up (see [5], [B] or [7]). Mechanical models 
of two planar systems, which have an unactuated cyclic variable, and all their shape 
variables are independently actuated, are considered in [8]. To deal with the stability 
of a walking gait for a biped, which is under actuated in single support, because it 
has point feet, in the papers [5], |f Oj and |f f ] . reference trajectories are defined for 
the actuated variables as functions of an undriven strictly monotone state variable. 
A complete characterization of all mechanical systems with underaction degree one is 
given in [12] . In [T3], a control law is proposed to stabilize the surge, sway and angular 
velocities of the hovercraft system. We can also note the very interesting thesis docu- 
ment [14] . which is devoted to nonlinear control, reduction, and classification of under 
actuated mechanical systems and in particular of high order under actuated systems. 
Then numerous mechanical devices associated with underactuation have been studied 
in literature. Furthermore this topic is far to be closed, because control design methods 
do not exist for many under actuated systems that are important for applications. 

This paper deals with the stabilization of two planar under actuated systems. The 
first system is the well-known straight beam-and-ball system. The ball is perfectly 
rolling without slide on the beam. Due to the complexity of this system, the stabiliza- 
tion and the tracking problem using a state or an output feedback have been considered 
by many researchers (see [15], [16], [17], [18] or [14]). In paper [15], tracking for this sys- 
tem was considered using approximate input-output linearization. Semiglobal stabiliza- 
tion of the straight beam-and-ball system using state feedback was addressed by [17] . 
In [16] . this system is stabilized using output feedback. The problem of global stabiliza- 
tion of the straight beam-and-ball system with friction was considered in paper [18] - 
The viscous friction is taken into account in our paper too. Semiglobal stabilization of 
this system, using fixed-point state feedback was addressed by |14j . 

The second system is an original circular beam-and-ball system. For each system, 
a control law, based on the linearized model and its Jordan form is designed. The 
saturation of the actuator is taken into account explicitly, so the control law is non- 
linear. This kind of control has been previously tested to stabilize a biped with point 
feet [T9], a one-link pendulum with flywheel [6], and to stabilize a two- link pendulum 
with flywheel [20] , The main difference between the straight beam-and-ball system and 
the circular beam-and-ball system is that the linear model of the second system has two 
eigenvalues in the right-half complex plane. Therefore, it is more difficult to stabilize 
the circular beam-and-ball system than the straight beam-and-ball system with only 
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one eigenvalue in the right-half complex plane. For the linearized model of the straight 
beam-and-ball system, the controllability domain, noted Q, and the basin of attraction, 
noted B, can coincide under a linear control law with restriction (see [2T], [5]). This 
property is not satisfied for the linearized model of the circular beam-and-ball system 
with two eigenvalues in the right-half complex plane. But the basin of attraction B can 
be made arbitrary close to the controllability domain Q (see [22]). 

We hope that our study is of theoretical interest and also has some pedagogical 
value. 

The paper is organized as follows: Section [2] is devoted to the straight beam-and- 
ball system. In Subsection 12.11 the equations of motion are written. The linear model 
is presented in Subsection 12.21 Subsections 12.31 12.41 12.51 are organized to introduce a 
control law with saturation, to get a basin of attraction B, which coincides with the 
controllability domain Q. The circular beam-and-ball system is studied in Section [3] 
In Subsection 13.11 the equations of motion are written. The linear model is presented 
in Subsection 13.21 Subsections 13.31 13.41 13.51 are organized to introduce a control law 
with saturation, to get a large basin of attraction B inside the controllability domain 
Q. Simulation results for the complete nonlinear systems are shown to illustrate the 
efficiency of the proposed control laws. Finally, Section [4] contains our conclusion and 
perspectives. 



2 Straight beam-and-ball system 

The straight beam-and-ball system consists of a straight beam and a ball on it, see 
Figure [T] The ball is rolling on the beam without slide. The point C\ is center of mass 
of the beam with its holder OA. The point C2 and and value r are center and radius 
of the ball. The point C2 is also the center of mass of the ball. 



2.1 Equations of motion 

Let mi and mi denote the mass of the beam with its holder OA and the mass of the 
ball, respectively. Let us introduce p\ and p2 the radii of inertia such that I\ = m\p\ 
and I2 = Tn,2p2 are respectively the inertia moment of the beam with its holder OA 
around the suspension point O and the inertia moment of the ball around its center 
C 2 ; let OCi = a and OA = I. 

Two generalized coordinates, the angular variables 6 and ip characterize the be- 
havior of this system. Position of the ball on the beam is defined also by the distance 
s = rip. Let F be the torque, which is directly proportional to the electrical current 
in the armature winding. By neglecting the armature inductance (in other words, the 
electromagnetic time constant in the rotor circuit), this torque can be written in the 
form (see [2"3"]): 

r = Cu u-c v e (1) 

where u is the voltage, supplied to the motor. The positive constants c u and c v for a 
given motor can be calculated by using the values for the starting torque, the nominal 
voltage, the nominal torque and the nominal angular velocity |23| . Product c v 6 is the 
torque of the back electromotive force. The torque of the viscous friction force in the 
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joint O (if it is taken into account) is also proportional to angular velocity 9. We will 
consider the following constraint, imposed on the voltage u: 

\u\ < urji u = const (2) 

The expressions for the kinetic energy K and the potential energy 77 are the fol- 
lowing (g is the gravity acceleration): 

2K = mipje 2 + m 2 [r V + (r + l) 2 ]d 2 + 2m 2 r(l + r)<p9 + m 2 (r 2 + p 2 2 )<Jp 2 

(3) 

77 = raigacosO + m 2 g\—rtpsinQ + (I + r)cos9] 

The equations of the mechanism motion can be derived, using Lagrange's method: 

[mipi + m 2 (r + I) 2 + m2r 2 tp 2 ] 8 + jn 2 r(r + l)tp + 2m 2 r 5 (pipO— 
— g[mia + m 2 (r + l)\sin9 — m2grtpcos9 = c u u — c v 6 

r(r + 1)9 + (r 2 + p 2 )</3 — r 2 ip9 2 — grsin9 = (5) 
If u — 0, system (|4]), ([5} has one unstable equilibrium state: 

9 = 0, ip = {s = 0), 6 = 0, ip = {s = 0) (6) 
2.2 Linearized Model 

Corresponding to the nonlinear equations ([4J , (J5J , the linear model of the motion near 
the unstable equilibrium state © is: 

[mip? + m 2 (r + Z) 2 J 9 + m 2 r(r + l)Cp - g[mia + rn 2 {r + 1)]9 - m 2 grip = c u u — c v 9 (7) 
r{r + 1)9 + (r 2 + p 2 2 )Cp - gr9 = (8) 



2.3 Kalman controllability 

The determinant of the controllability matrix (see [24]) for the linear model (0, (|8} is 
not null, if and only if: 

r 2 g 2 [{2r 2 + pi) pi +r i ] (9) 

Thus, inequality (|9j) is valid, if r ^ 0. If r = 0, then the ball becomes a material point 
and we do not consider this case. Thus, the linear model of the straight beam-and-ball 
system is always controllable. 
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2.4 Spectrum of Linear System 

The state form of system l[7 |) . (8 ) . using the state vector x — (9, ip, 8, ip) 1 



x = Ax + bu 



2 > 



D~ E D~ 



-ov 




x + 



D 






-1 ( Ou 





(10) 



The notations 02x2 and ^2x2 define a zero matrix and an identity matrix, respectively. 
The expressions of matrices D and E are 



D 



m\p\ + m 2 {r + I) 2 m 2 r(r + I) 
r(r + l) r 2 + pl 



(11) 



E = g 



mio + m,2{r + I) rri2T 
r 



Introducing a nondegenerate linear transformation x — Sy with a constant matrix 
S, it is possible to get the well-known Jordan form of the matrix equation (|10|l 



y = Ay + du 



(12) 



where 



A = S~ L AS 



Ai \ 

V A 4 / 



d^S- 1 b=[d l ] T (< = !,.. .,4). (13) 



Here, Ai, A4 are the eigenvalues of the matrix A. They are the roots of the charac- 
teristic equation of system (O, (|8j): 



432 
aoA +aiA + a 2 A + a^X + 04 = 



(14) 



with 



a = detD > 0, ai = c^(r 2 + p 2 ) > 0, a2 = ms(r + l)(r 2 - p\) - ra\a{r 2 + p 2 ), 
a 3 = 0, a 4 = detE = -m 2 g 2 r 2 < 0. 



If all physical parameters of the studied system are known, matrix S of the transfor- 
mation x = Sy can be calculated. 

According to the theorem of Routh-Hurwitz (see [25]), equation (|14[) has one root 
in the right-half complex plane and three roots in the left-half complex plane (see 
also [14.). This assertion does not depend on the sign of the coefficient a 2 . Of course, 
the unique root in the right-half complex plane is located on the real axis. 

In Section |3l we consider the ball on the circular beam and use linearized model 
in the same matrix form (|10[1 as for the ball on the straight beam, but with different 
submatrices D and E. 
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2.5 Problem Statement 

Let x = (here is a (4 x 1) zero-column) be the desired equilibrium state of system 
(|10|l . Let us design the feedback control u(x) to stabilize this equilibrium state x — 0, 
under constraint @. In other words, we want to design an admissible (satisfying the 
inequality ([2])) feedback control |w(:c)| < uq to ensure the asymptotic stability of the 
desired state x — 0. Let W be the set of piecewise continuous functions of time u(t), 
satisfying inequality @. Let Q be the set of the initial states x(0) of system {TO}, from 
which origin x = can be reached, using admissible control functions of time u(t). 
In other words, system (|10[) can reach the origin x = with the control u(t) G W, 
only starting from the initial states x(0) € Q. Set Q is called controllability domain. 
If the matrix A has eigenvalues with positive real parts and the control variable u is 
restricted, then the controllability domain Q for system (|lUp is an open subset of the 
phase space X (see [21], [6]). 

For any admissible feedback control u = u(x) with saturation |«(a;)| < uq the 
corresponding basin of attraction belongs to the controllability domain: B C Q. Here, 
as usual, B is the set of initial states x(0), from which system (|10|l . with feedback 
u — u(x) asymptotically tends to the origin point x = as t — > oo. 

In the following section, a control law will be presented for the straight beam-and- 
ball system to get a basin of attraction B, which coincides with the controllability 
domain Q: B = Q . 

2.6 Feedback Control for the straight beam-and-ball system 

A control law is proposed here to stabilize the straight beam-and-ball system with 
basin of attraction as large as possible. 

2.6.1 Control design 

Let Ai be the real positive eigenvalue, Re\ < (i — 2, 3, 4) and let us consider the 
first scalar differential equation of system H12[) corresponding to eigenvalue Ai, 



System (jTUJ) , is a Kalman controllable system, therefore scalar d\ 7^ 0. The controlla- 
bility domain Q of the equation (| 15p and consequently of system ()12[) is described by 
the following inequality (see [21], [6]) 



The instability of the coordinate y\ can be "suppressed" by a linear feedback con- 
trol, 



Vi = Aiyi + diu 



(15) 




(16) 



u = 72/1 



(17) 



with the following condition, 



Ai +di7 < 



(18) 



For system (| 10p under the feedback control ()17p with inequality (|18|) . only the pole 
Ai is replaced by a negative pole Ai + ^17. The poles A2, A3, A4 do not change. 
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If constraint ([2| is taken into account, the linear feedback control (|17[) becomes 
with saturation, 

WQi if 72/1 > u 

= = s jyi, if \ivi\ < «o (19) 

-uo, if 72/1 < -WO 
The unit of coefficient 7 is volt. 

It is possible to see that if \yi\ < di|tto/Ai, then under condition (|18|l the right 
part of equation (|15p with the nonlinear control (|19p is negative when y\ > and 
positive when y\ < 0. Consequently, if |j/i(0)| < |di|uo/Ai, then the solution yi(t) 
of system (|15[) . ()19|) tends to as t — > 00. But if yi(t) —¥ 0, therefore, according to 
expression ()19p . w(t) — ► as t — > 00. Therefore, the solutions yi(t) (i = 2,3,4) of 
the second, third and fourth equations of system (|12[) with any initial conditions 2/j(0) 
(z = 2, 3, 4) converge to zero as t — > 00, because ReXi < for i = 2, 3, 4. Thus, under 
the nonlinear control H19[) and with inequality (|18|l . the basin of attraction _B coincides 
with the controllability domain Q (see [2T], [5]): B = Q. So, the basin of attraction B 
for system (|10l) . (|19[) is as large as possible and it is described by inequality (|16l) . 

Note that the variable yi depends on the original variables from the vector x, 
according to the transformation x = Sy or y — S~ 1 x. Due to this, formula (|19p defines 
the control feedback, which depends on the vector x of the original variables. If the 
matrix S is calculated, then all coefficients of the designed control can be defined. Only 
the constant 7 is an arbitrary multiplier, but it has to satisfy inequality (1181) 

Thus, linearizing nonlinear system @, (JSJ), (|19p near the equilibrium state we obtain 
a system, which is asymptotically stable. Using Lyapounov's theorem (see [26]), we 
conclude that equilibrium (|6| of the nonlinear system ((4J), (J5j) is asymptotically stable 
under control ()19p with some basin of attraction. In the next Subsection, numerically 
we find the upper bounds of the initial values of some variables, which can be handled 
for the linear and nonlinear models under the designed control. 



2.6.2 Numerical results 



Let 

mi = 1.0 kg, m,2 — 0.2 kg, g = 9.81 m/s 2 , 
r = 0.05 m, I = 0.2 m, a = 0.15 m, Pl = 0.2179 m, p 2 = 0.1414 m, (20) 

c u = 0.007 N-m/V, c v = 0.0001 N-m/s, u = 19 V. 

In open- loop the poles of the linear system (| 10p (the roots of equation 1)141) ) with 
parameters (|20p are: 

Ai = 5.7202, A 2 = -5.7218, A 3 , A 4 = -2.8-10~ 7 ± 1.0558i, (21) 

Now we can use inequality (1161) to evaluate the basin of attraction B for system 
(flU)) , (TT91 . If 8(0) = 9(0) = ip(0) = 0, the upper bound of the initial angles ip, which 
can be handled for the linear model (| 10p is ip(0) = 77.679°. The corresponding initial 
distance s(0) is equal to 0.0678 m. This value for the distance s is close to the value 



8(0) = ^ 

m 2 g 



(22) 
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With 9 — product sm 2 t/ is the torque about joint O of the gravity force of the ball 
(see the nonlinear equations (|4]), (|5]) and the linear equations Q, p|). the product 
c u uq is the torque (maximal as possible) developed by the motor in static. Thus, the 
point 

9 = 6 = a = 0, a = ^ (23) 

is the equilibrium state (unstable) for our system (nonlinear @, © and linear l[7]l. 
p)l). It is easily to see that the equilibrium point (|23[) is located on the boundary of 
the controllability region H16[) . Simulation shows that, if 

6(0) = 0(0) = s(0) = 0, s(0) > (24) 

then it is not possible to bring the nonlinear system @, (jSJ under control (|19l) to the 
equilibrium (|6|; but it is possible to do that, if s(0) < . Furthermore, we think 
there is no an admissible control |w(x)| < mo to bring system (|4]), ([5} to the equilibrium 
point (0 from the initial states (|24p . This opinion is based on the numerical studies 
and physical feeling. We do not prove here corresponding assertion strictly. 

The eigenvalues A3, A4 are very close to the imaginary axis (see (|21) l1 and therefore 
under the control (|19|l . the transient process is very long. Let us take into account a 
viscous friction in the joint O defined by the torque fd. The consideration of the torque 
fd of the friction force is equivalent to the consideration in equation JT]) of the term 
(c v + f)6 instead of the term c v 9. With / = 0.4 N-m-s for example the poles of the 
corresponding linear system (|10|l in open-loop are: 

Ai = 3.4001, A 2 = -10.0181, A3, A 4 = -0.1041 ± 1.0297i (25) 

The technique of the feedback control design with a viscous friction (with new poles 
H25[) ~) remains the same exactly. And the structure of this control remains the same - 
(|19l) . Under the control law (| ltjp with new coefficients, the transient process converges 
to the equilibrium state ([B]) faster than without friction. Using inequality (|16|) . or the 
equality (|22l) we get of course the same value s(0) as above without friction. So, we can 
use formula ()22|l for the linear and nonlinear systems to calculate the upper bound of 
the initial distances s, which are possible to stabilize the equilibrium state ©. 

Figures [5] and show a numerical test with an initial tilt ip(0) = 77.65° for the 
nonlinear system (|4]), ((5]) with the coefficient / = 0.4 N-m-s under the control law (11911 
with 7 = —122. The voltage, supplied to the motor, is shown in Figure [3] The limit 
value uq — — 19 V is reached at initial time. 

Let F be the reaction force, applied to the ball orthogonally to the beam in their 
contact point. The following formula for this force holds: 

F = m 2 \g cos9 - (I + r)6 2 - 2np6 - rtp§] (26) 

If the reaction force F becomes negative, then the ball loses contact with the beam 
and our model (with contact) becomes false to describe the physical process. In the 
numerical experiment, presented in Figures [2] and [3] the force F is always positive. 
This force is shown in Figure 0] 

If 95(0) = <£(0) = 9(0) — 0, then, using inequality H16[) . the upper bound of the 
initial tilts of the beam, which can be handled, for the linear model (j 10|) with the 
friction is 9(0) = 3.61°. The computations show that the upper bound of the initial 
tilts for the nonlinear system @, (© under control (|19p is 6(0) — 3.64°. So, this value 
is little more important than for the linear system (|10[) under the same control (|19[) . 
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3 Circular beam-and-ball system 

The circular beam-and-ball system consists of a circular beam with the center C and 
the radius R and a ball on it with the center C 2 and the radius r, see Figure [5] The 
point C\ is the center of mass of the beam with its holder OA. 



3.1 Equations of motion 

Here the same notations are used, that for the straight beam-and-ball system. 

Let mi and m 2 denote the mass of the beam with its holder OA and the mass of 
the ball, respectively. Let p\ and p 2 be the radii of inertia respectively of the beam 
with its holder OA and of the ball; let OG\ = a and OA — I be. 

The generalized coordinates are the joint variable 8 and the angle variable ip. Po- 
sition of the ball on the beam is also defined by distance s — rip. The relation between 
the angle ip and angle tp is: 

rip = Rtj) 

Let us assume that the motor is the same that for the straight beam-and-ball 
system with torque ^ and constraint (J2]). The constants c u , c v and ito are the same. 

The expressions for the kinetic energy K and the potential energy 77 are the fol- 
lowing: 



2K= {m lP \+m 2 [{R + r) 2 + {I - R) 2 + 2(R + r)(l - R)cos^]} 2 + 
+m 2 [(R + r) 2 + ^l]^ + 
+2m 2 [(R + r) 2 + (R + r)(l - R)cos^] 9^ ( 27 ) 

77 = [mi a + mzil — R)]gcos6 + m 2 g(R + r)cos(fytp + 8) 

The equations of the mechanism motion are derived, using Lagrange's method: 

[mip\ + m 2 {r 2 + I 2 + 2rlcos^-) + 2m 2 R(R + r-l)(l- cos^)] 8+ 
+m 2 r(l + £) [7? + r + (I - R)cos^] (p + m 2 r(l + - l)(28 + %)<prinZ$- 

— g[m\a + m^il — R)]sin6 — m 2 g(R + r)sin(8 + ^) = c u u — c v 8 



r(l + %) [R + r+(l- R)cos^} 8 + [p 2 +r 2 (l + £) 2 ] Cp+ 
+(1 + fc)(J - R)8 2 sin^ - gr(l + £)«n(0 + = 

If u — 0, system (|28p . (|29|l has one unstable equilibrium state (|6]). 



(28) 



(29) 
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3.2 Linearized Model 

Linearizing the equations (|28l) . (1291) near the unstable equilibrium state IJBJ, we get the 
following model: 

[rn-ipl + m 2 (r + I) 2 ] 6 + +m 2 r (l + £) (r + 

(30) 

— g[mia + m 2 (r + Z)]0 - ra 2 #(r + -R)^ = c n u - c«0 
r (l+ ^) {r + i)e + [p 2 +r 2 (l + £) 2 ] 0- 5r (l+ + = (31) 

3.3 Kalman controllability 

The determinant of the controllability matrix for the model (|3U|) . (|31|l is not null, if 
and only if: 

Rr 2 (R-l) + R 2 pl+r 3 (R-l)^0 (32) 
If r = 0, then the ball becomes a material point and p 2 = 0. In this case, instead 

of inequality (|32p the equality is correct. However, we do not consider a material point 

on the beam and therefore assume r ^ 0. 

Let r / 0, but the mass of the ball is concentrated in its center (p 2 = 0) and the 

suspension point O coincides with the curvature center C of the circular beam (R = I). 

In this case, inequality (|32[) is not satisfied and the linear system is not controllable. 

Consider the controllability of the original nonlinear system (|28|) . (|29l) in the case p2 = 

and R = I. Introduce the angle a = 8 + The nonlinear system (|28l) , (I29p becomes: 

m\p 2 Q — rrtigasinO = c u u — c v 9 (33) 

(R + r)a — gsina — (34) 
The equations (|33[) and (|34[) are separated. The control u has no action on the angle 
a and system (J33J) , (|34f) is not controllable. 
Inequality (|32p is satisfied, if 



R 2 2 
(-ft + r)r z 



and we will consider only this case. 



3.4 Spectrum of Linear System 

The state form of system (|30[1 . (|31[) can be presented in the same matrix form (110[) as 
for the straight beam, but with the following submatrices D and E: 



D 



'm 1 pl + m 2 {r + l) 2 m 2 r(l + £)(r + Z) ' 
. r(l + ^)(r + I) p| + r 2 (l + £) 2 , 



(36) 

'mia+m 2 (r + /) m 2 r(l + -^)\ 
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Introducing a nondegenerate linear transformation x = Sy with a constant matrix 
S, we can get the Jordan form similar to (|12|) . (|13|l . 

The characteristic equation of system (130 j> . (|31j> has form (|14|) with 

a = detD>0, a x = c v \p\ + j^{R + r) 2 j > 0, 

a 2 = -mis [(/al + r 2 )aR 2 + (2ar + p\)Rr 2 + (ar + pf)r 3 ] - 



Inequality (|37[1 is satisfied, if the radius R of the circular beam is sufficiently small (the 
curvature of the beam is sufficiently large). But we have not to forget condition (|32|l 
(or ((35j ) °f controllability. 

Under condition (|3Tf) . the coefficient 04 is positive. Using the theorem of Routh- 
Hurwitz (see [25]), we can conclude that the characteristic equation (j 14f) has two roots 
in the right-half complex plane and two roots in the left-half complex plane. This 
conclusion does not depend on the sign of the coefficient a 2 . 

3.5 Problem Statement 

We will consider the same problem, as before for the straight beam-and-ball system. 
We want to design an admissible (satisfying inequality p|) feedback control to ensure 
the asymptotic stability of the state x = with a large basin of attraction for this 
equilibrium state. 

3.6 Feedback control for the circular beam-and-ball system 

A feedback control law u(x), satisfying inequality {2), is proposed here to stabilize the 
circular beam-and-ball system with a large basin of attraction. Under condition (|37[) . 
the linear model of the system has two eigenvalues Ai, A2 in the right-half complex 
plane and two eigenvalues A3, A4 in the left-half complex plane. 

3.6.1 Control design 

Let Ai and A2 be the real positive eigenvalues, and let us consider the first two scalar 
differential equations of system (|12[1 . (|13p for the circular beam, corresponding to these 
eigenvalues Ai and A2: 



-m 2 g [{r + l)R 2 (r 2 

2 

a 3 = -c v gjp(R + r) < 0, 




We assume that 



rri\a + 1x12(1 — R) > 



(37) 



Vl = + diu, y 2 = A 2 j/2 + d 2 u 



(38) 



Under condition (|35[) system (| lUp for the circular beam is Kalman controllable. 
Therefore, subsystem (|38|) is controllable too (see 24 ) and d\ 7^ 0, d 2 7^ 0. The 
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controllability domain Q of the equations (|38|l , and consequently of system (|12|l , is an 
open bounded set with the following boundaries (see |27p 



W ( T ) = ±2£ 

(39) 

y2 (r) = ±*^a (2e~ X2T -1) (0 < r < oo) 

If the system has two complex poles in the right-half complex plane, then instead 
of (|39[) we will get other formulas (see \21\). 

Set Q belongs to the rectangle, defined by inequalities: 

\vt\<\di\^, \y 2 \<\d 2 \% 
The boundary of the controllability region Q has two corner points (see Figure [6]): 

Vi = -di%, v» = -<h%; 

(40) 

Vi=d!^, y 2 = d 2 ^ 

These points (|4U|I are the equilibrium points of system (|38p under the constant controls: 

u = ±uq (41) 

We can "suppress" the instability of the state yi — 0, y 2 — by a linear feedback 
control, 

« = fciyi + k 2 y 2 (42) 

with fci = const and A;2 = const. It is shown in paper [22] that using a linear feedback 
(142|l with saturation (7 = const): 

uq, if j(kiy\ + k 2 y 2 ) > u 

l{k\yi + k 2 y 2 ), if \j(kiyi + k 2 y 2 )\ < u (43) 

-wo, if l{kiyi + k 2 y 2 ) < -uq 

the basin of attraction B can be made arbitrary close to the controllability domain Q. 
The straight line crossing two points (|40|l is the following: 

kiyi + k 2 y 2 = 

with 

fci = -£, fca = f (44) 

A2 Ai 

If 

sign'y = sipn [did2 (^1 — ^2)] (45) 



and I7I — y 00, then the basin of attraction B of system (|38[) under the nonlinear control 
(143 II with coefficients (|44[) tends to the controllability region Q. Consequently, using the 
coefficients (|44[) . the basin B can be made arbitrary close to the domain Q. If I7I — 5* 00, 
control (|43fl tends to the bang-bang control. 



13 



The solutions yi(t) and y^it) of system (|38|) . (|43p tend to as t — v oo for the initial 
values 2/i(0), 2/2(0), belonging to the basin of attraction of system (138[) . (|43|l . But if 
2/i (t) — > and 2/2 (i) - ► 0, then, according to the expression (|43[1 . u(t) -> as t 00. 
Therefore, solutions ys(t), yn(t) of the third and fourth equations of system (|12p with 
any initial conditions 2/3(0), 2/4(0) converge to zero as t -> cxd, because ReX^ < 0, 
Re\4 < 0. Thus, under control (|43p with coefficients (|44fl . the basin of attraction 
of system (|12|l . (|43p is described by the same relations, which describe the basin of 
attraction of system 1)380. (143 1) . 

The variables 2/1 and 2/2 depend on the original variables from the vector x, ac- 
cording to the transformation y = S~ 1 x. Due to this, formula (143 } defines a nonlinear 
feedback control, which depends on the vector x of the original variables. If the matrix 
S is calculated, then all coefficients of the designed control can be found. Only the 
constant multiplier 7 is an arbitrary one; but it has to satisfy relation (|45p and to be 
sufficiently large in modulus. 

Thus, linearizing the nonlinear system (|28p . (|29p . (I43p near the equilibrium state we 
get the system, which is asymptotically stable. Then according to Lyapounov's theorem 
(see 26 ), the equilibrium state ((6} of the nonlinear system (128|) . (I29p is asymptotically 
stable under control (143 1 with some basin of attraction. In the next Subsection, nu- 
merically we find the upper bounds of the initial values of some variables, which can 
be handled under designed control. 

3.6.2 Numerical results 
Let 

mi = 1.0 kg, rri2= 0.2 kg, g = 9.81 m/s 2 , 

(46) 

r = 0.05m, R = 0.8m, I = 0.2m, a = 0.15m, p x = 0.2646m, p 2 = 0.1414m 

In open-loop the poles of the linear system (|10p with parameters (|46p are: 

Ai = 4.89589, A 2 = 0.46516, A 3 = -4.89706, A 4 = -0.46523 (47) 

Using formulas ((39}, the controllability domain Q for system (|38l) is designed. It is 
bounded in Figure [ll] by dashed line. Its boundary contains the corner points (|40l) . 
Using the linear model (|30p , (|31fl , we can define the following equilibrium points under 
controls (|41[l with the original variables: 

1 g[mia+mM— R)\ ' 

(48) 

= (s = -R6), <p = (s = 0). 

Points (|48D are located on the boundary of the controllability region. 

Using the nonlinear model (1281) . (|29p . instead of (|48p we get the following expres- 
sions: 

6 = Tarcsin g[mia c + ^° {l _ R)l , 6 = 0, 

(49) 

<p = -§9 (s = -R6), <y5 = (s = 0). 

In equilibriums (|49fl . the ball is located on the highest point of the circular beam and 
consequently in this point of contact between the ball and the beam the tangent to the 
beam is horizontal. Remind, if u = 0, we have the equilibrium state ([6]). 
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The basin of attraction B for system (138p under the control (143 1 is shown in the 
same Figure [6] Its boundary is drawn in Figure [6] by solid line. This boundary is 
the periodical motion (cycle) of system (|38p , (|43p . This cycle is computed, using the 
backward motion of system (|38p , (|43)1 from a state close to the origin yi — yi — 0. 
The basin B depends on the coefficient 7. We show in Figure [6] the basin of attraction 
B, with 7 = 4000. If the coefficient 7 is smaller, then the basin of attraction is smaller 
too. 

In simulation, the control law (|43p is applied to the nonlinear model (|28|1 , (|29l) . 
Figure [7] shows the graphs of the angular variables and ip. These graphs are designed 
for the initial angle <p(0) = 70.39° and 0(0) = 0(0) = <p(0) = 0. This value <p(0) = 
70.39° is close to the upper bound of the initial angles <p(0), which are possible to 
stabilize the equilibrium state ([SJ. The corresponding initial distance s(0) is equal to 
0.061 m. No oscillations appear during the transient process, because matrix A does 
not have complex poles. The voltage, supplied to the motor, is shown in Figure [8] The 
limit value Uq = —19 V is reached at initial time. No oscillations also appear in the 
graph of voltage u(t). 

The following formula holds for the reaction force F, applied to the ball orthogonally 
to the beam in their contact point: 

1 + ^) - (R + r) (9 + ^ j + (R — l)9 2 cos T ^- -(R- l)9sin^ 

(50) 

In the numerical experiment, presented in Figures and [51 force F is always positive. 
This force is shown in Figure [9] 

Calculating the values and s with first and third formulas in (|49fl . we obtain 



F — mi 



g cos 



9 = 0.469, s = -0.375 m (51) 

Consider the initial velocities 0(0) = 0, s(0) = and let s(0) = — R9(0) be (see 
third equality in (|49p ). Simulating the nonlinear system (I28p . (I29p under control (I43f) 
(7 = 4000), we get values which are close to the boundary of the attraction basin 

0(0) = 0.397, s(0) = -0.318 m (52) 

These values (|52f) equal to values ()5ip divided by 1.18. Under the nonlinear control law 
Q43p with 7 = 8000 we come to the values 

0(0) = 0.430, s(0) = -0.344 m 

These initial conditions equal to values Q5ip divided by 1.09 and they are closer to 1)5 H 
than values ()52p . Our numerical experiments show that possible for stabilization initial 
values 0(0), s(0) tend to values ()5ip as 7 — > 00. Thus, formulas (|49p can be used to 
evaluate the basin of attraction for the original nonlinear system ()28p . ()29l) . 

We think there is no admissible control |w(a;)| < uq to bring system 1)28)1 . ()29|) to 
equilibrium © from the initial states 

0(O) = ^(O) = O, |0(O)|>armn g[min ^° (i _ fl)1 , 
V (0) = -#0(0) («(0) = -R9(0)). 



This hypothesis is similar to the corresponding hypothesis for the straight beam-and- 
ball system (see subsection 12.6.21) . 
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4 Conclusion 



In this article, we consider the well known straight beam-and-ball system and an orig- 
inal circular beam-and-ball system. The problem of stabilization of unstable equilibri- 
ums of these systems is studied. The model linearized near the unstable equilibrium 
of the straight beam-and-ball system has one unstable mode. The difficulty is greater 
to stabilize the circular beam-and-ball system, because its linear model has two un- 
stable modes. For each system we use the Jordan form of the linear model to extract 
the unstable part and to stabilize the equilibrium. Considering the restriction on the 
voltage of the motor, the objective is to get a large basin of attraction. The designed 
feedback control contains the unstable Jordan variables only. All parameters of this 
control are defined up to a constant multiplier. Simulation results for the complete 
nonlinear systems are shown. These results are close for linear and nonlinear systems. 
All the numerical results, obtained in this paper for nonlinear systems, are realistic 
and illustrate the efficiency of the designed control laws. Using the described above 
approach for both cases it is easily to take into account the friction forces also. Testbed 
devices can be now imagined to test the designed control laws experimentally. The 
original circular beam-and-ball system will be interesting for demonstrations, devoted 
to education and to investigate new nonlinear control laws. 



5 List of Captions 



1. Figure [T] Diagram of the straight beam-and-ball system 

2. Figure [2] Stabilization of the Straight beam-and-ball system, 8 — > and <p — > (in 
radians) . 

3. Figure Voltage u(t), supplied to the motor for the stabilization of the Straight 
beam-and-ball system. 

4. Figure [4] The reaction force F(t), applied to the ball during the stabilization pro- 
cess. 

5. Figure [S] Diagram of the circular beam-and-ball system. 

6. Figure [6] Controllability domain Q (dashed line) for system (|38[) and basin of 
attraction B (solid line) with 7 = 4000. 

7. Figure [7] Stabilization of the Circular beam-and-ball system, 9 — ¥ and ip — s> (in 
radians) . 

8. Figure [8] Voltage u(t), supplied to the motor for the stabilization of the circular 
beam-and-ball system. 

9. Figure [§1 The reaction force F(t), applied to the ball during the stabilization pro- 
cess. 
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Fig. 1 Diagram of the straight beam-and-ball system. 
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Fig. 2 Stabilization of the Straight beam-and-ball system, 6 — » and <p — > (in radians). 
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Fig. 5 Diagram of the circular beam-and-ball system. 
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Fig. 7 Stabilization of the Circular beam-arid- ball system, 6 — > and <p — ¥ (in radians). 
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Fig. 8 Voltage u(t), supplied to the motor for the stabilization of the Circular beam-and-ball 
system. 
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Fig. 9 The reaction force F(t), applied to the ball during the stabilization process. 
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